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ABSTRACT 



We present line lists for ammonia (NH3) at high temperatures obtained by 
recording Fourier transform infrared emission spectra. Calibrated line lists are 
presented for twelve temperatures (300 - 1300°C in 100°C intervals and 1370°C) 
and each line list covers the 740 - 2100 cm~^ range, which includes the majority 
of the U2 umbrella bending mode region (11 /xm) and the 1/4 asymmetric bending 
mode region (6.2 jim). We also demonstrate the useful technique of obtaining 
empirical lower state energies from a set of spectra recorded at different sample 
temperatures. Using our NH3 spectra, we have estimated lower state energies 
{Elow in cm~^) and our values have been incorporated into the line lists along 
with calibrated wavenumbers (z> in cm~^) and calibrated line intensities {S' in 
cm/molecule) for each line. We expect our hot NH3 line hsts to find direct 
application in the modelling of planetary atmospheres and brown dwarfs. 

Subject headings: Emission spectra; Molecules(ammonia); Stars(cool, low mass); 
Brown dwarfs; Atmospheric models; Line lists 
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INTRODUCTION 



High-resolution infrared laboratory s pectra are an essenti al element in nu merous fields 



of sci ence such as atmospheric chemistry (IBernath et al 



20051 ) and astronomy (ISwain et al. 



20081 ). These high- resolution spectra allow molecular emission and absorption features to 



be accurately identi fied and assigned; for example the detection of 'Water on the Sun' 



( 1 Wallace et al. 



19951 ) was made possible by the use of high-resolution hot water (H2O) line 
lists derived from experimental spectra. The H2O molecule is able to exist in sunspots 
(typically 3,000 - 4,500 K), and observations showed that many previously unassigned lines 
were in fact due to the ubiquitous H2O molecule. 



Molecules form in 'cool' sources ( iBernathl l2009l ) and progress in infrared astronom y 



has led to t 



le ide ntification of nu merous molecules in cool stars (IWing fc Ford 



AUard et al. 



19971 ). bro wn dwarfs ( 



Gebal 



recently, in exoplanets (ISwain et al. 



2009 



e et al. 



2002 



Sharpe fc Burrows 



Seager fc Deming 



20101 ). 



1969; 



20071 ) and, most 



After H and He, C, N and O are the most abundant elements, so in cool objects one 
anticipates that methane (CH4), NH3, and H2O will be abundant. Although the spectra of 
the fundamental infrared bands of 'cold' (i.e. room temperature) CH4, NH3, and H2O have 
been understood for many years, the situation for overtones, hot bands and combination 
bands is not so sanguine. For H2O it is only recently that high resolution experimental 
and theoretical line lists suitable for computing the spectral energy distributions (SEDs) of 
astronomical obiect s with temperatures in the 300 - 4000 K range have become available 



( ICoheur et al. 



20051 : 



Barber et al 



20061 ). For 'hot' CH4, there are only some rather sparse 



experimental observations and no satisfactory quantitative the oretical understanding. 



For hot NHcj there has been some rec ent theoretical progress ( lYurchenko et al. 



2011 



Zobov et al. 



2011 



Huang et al. 



2011al Jbl) and we report here on an extensive set of 



experimental observations in the thermal infrared. 
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The NH3 molecule is one of the simplest polyat omic mo 



ecules with four atoms forming 



20051). NH3 has 3A^ - 6 = 6 



a trigonal pyramidal structure with C^y symmetry (iBernathI 
fundamental vibrational modes, of which two are doubly degenerate. There are therefore 
four fundamental vibrational frequencies i^i(cii) at 3336.2 cm~^ (N-H stretch), 1^2(01) at 
932.5 cm~^ (umbrella mode), i^sie) at 3443.6 cm^^ (N-H stretch) and z/4(e) at 1626.1 cm~^ 
(bend), of which z/3 and z/4 are degenerate (Table [1]). There is an additional complication in 
the spectroscopy of NH3 due to the low barrier to 'inversion'. Although NH3 is non-planar 
at equilibrium the three hydrogen atoms can flip through the plane containing the N atom 
and perpendicular to the C3 symmetry axis. 
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Table 1. The vibrational modes and rotational constants of NH3 ( lHerzberg]|l99ll ). 



Mode Symmetry Frequency (cm ^] 



Type 



Rotational Constants 



Z/4 



ai 3336.2 (3337.2* 

Qi 932.5 (968.3* 

e 3443.6 (3443.9* 

e 1626.1 (1627.4* 



Symmetric Stretch 
Symmetric Bend (Umbrella) 
Asymmetric Stretch 
Asymmetric Bend 



Bo 
Co 
ro 
a 



9.9443 cm-i 
6.196 cm~^ 
1.0173 A 
107.8 ° 



*Molecular inversion causes the modes to be doubled. 
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If the H atoms were labeled this would correspond to a right-handed form of the 
molecule converting into a left-handed form. Hence, there are two forms (frameworks) of 
the NH3 molecule (energy levels are conventionally labeled as s and a or + and -) which 
can interconvert. Indeed, the interconversion frequency at 24 GHz {s — a transition for 
J = 1, K = 1) can be u sed by radio astronomers to monitor NH3 in dark interstellar clouds 



teenson fc Myers 



1989I ). In general, this inversion motion doubles vibration-rotation lines 



(each of the two interconverting forms has slightly different energy levels) and leads to the 
two slightly different band origins listed in Table [1] 

NH3 is an oblate symmetric top {B = 9.9443 cm^^, C = 6.196 cm~^) and the 
fundamental modes are either parallel {ui and z/2, ai — ai) or perpendicular (z/3 and u^, 
e — ai) depending on whether the infrared transition dipole moment is oscillating parallel 
or perpendicular to the C3 symmetry axis. The rotational quantum numbers J and K, 
describe the quantization of the total angular momentum (excluding nuclear spin) and 
its projection onto the symmetry axis, respectively. Rotational selection rules {AK = 0, 
A J = 0, ±1) for a parallel transition lead to a simple P, Q, R structure (A J = —1,0, 1) 
with a small K splitting. The selection rules for perpendicular transitions (Ai^ = 0, ±1, 
A J = 0, ±1) lead to a more complex pattern with widely spaced K struct ure, organized 



into K sub-bands each with different K' — K" values and P, Q, R structure (iBernath 



20051 1. 



Over the y ears, numerous studies have been undertaken to assign the comp 



Cottaz et al 



NH3 spectrum (Povnter fc Margolis 



200G 



Lees et al. 



1983 



1984 



Sasada et al. 



1992 



Kleiner et al 



icated 



1999 



20081 ). Many studies are of cold NH3 and are useful when 



studying planetary atmospheres in our solar system (IBrown fc Margolis 



has been done in determining empirical lower state energies 



well as on the rotation-inversion spectra fiSasada et al. 



1992 : 



1 99 61) . Some work 



Brown fc Margolis 



Yu et al 



19961 ) as 



2OI0I ). Recently there 



has been considerable progress in the variational calculation of vibration-rotation energy 
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levels an d transitions from adiusted a b initio potential energy surfaces (IHuang et al.ll2008l . 



2011a 



b 



Yurchenko et al 



2009 



at 300 K (j Yurchenko et al 



2011 



). Calcula ted line lists suitable 



20091 ) and 1500 K flYurchenko et al. 



'or spectra of ammonia 



201 ll ) are now available. 



Cool objects are less luminous than hot stars and emit most of their radiation in the 
infrared. High resolution infrared spectra of these objects are best recorded with large 
10 m class telescopes. In 1995 the fi rst brown dwarf, Gliese 2 29B was identified by the 



presence of methane in its spectrum (jOppenheimer et al 



19951 ). Obj ects with masses les s 



than 0.0075 solar masses are unable to fuse hydrogen in their cores (IBurrows et al. 



20011) 



and are called brown dwarfs. Brown dwarf atmos pheres are cool en ough (typically 700 - 



2,000 K) to sustain a rich molecular environment ( iKirkpatrickl 120051 ) . Since the discovery 



of Gliese 2 29B, the L and T c lasses have been introduced to describe objects cooler than 



M- dwarfs (jLeggett et al 



metal hydrides such as FeH (IHargreaves et al 



20011). L-dwarfs display ii ear infrare d electronic transitions of 



2010) and CrH flKirkpatrick et al 



have weak TiO and VO transitions which are characteristic of M-dwarfs (IBochanski et al 



1999[) and 



20071 ). In contr ast, T-dwarfs like Gli ese 229B have prominent overtone transitions of hot 



Telescope, 



H2O and CH/i (IBurgasser et al. 



Gushing et al. 



20061 ). Based on observations taken from the Spitzer Space 



( I2OO6I ) identified prominent H2O, GII4 and NH3 absorption bands 



between 6-14 fim, including the NH3 1/2 band near 11 fim. It is predicte d that a new coole r 



20051) 



classification, a Y-dwarf (< 700 K) is required below the T-dwarf class (IKirkpatrickl 
with the p rediction that strorig NH.s absorptions w ill help to differentiate the Y-dwarfs from 



T-dwarfs (jLeggett et al. 



2007 : 



Delorme et al. 



20081). 



Numerous unassigned features are present in the observed SEDs of brown dwarfs. 
Due to a lack of experimental data recorded in the temperature ranges appropriate for 
these objects (or reliable theoretical predictions), the best option is to e xtrapolate roora 



temperature 'cold' line parameters such as the HITRAN 2008 database (iRothman et al 
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20091 ) to the temperature of interest (typically 1000 K for sub-stellar brown dwarfs). For 
this reason, these extrapolated spectra rarely match observations satisfactorily (see Figure 
H]) due to the lack of hot bands and highly-excited rotational transitions which contribute 
significantly to the SED in these temperature ranges. 

Along with brown dwarfs, it has only recently become possible to detect planets 
ou tside of the Solar Syste m (exoplanets). 51 Pegasi b was the pioneering discovery 



by 



Mayor &: Quelozl (119951 ) that paved the way for many more exoplanet detections. 



with over 500 discovered to date ( |http:/ /exoplanet. eu/ ). The vast majority of these 
exoplanets are large gas giants close to their parent star with high temperature gaseous 
atmospheres t ypically referred to as 'ho t Jupiters'. The new technique of 'transit 



spectroscopy' ( 



Char 



jonneau et al. 



(ICharbonneau et al 



20001 ) has been used to probe exoplanet atmospheres 



2OO2I ) and familiar molecules such as H2O, CH4, CO and CO^ have all 



been shown to be p resent in examples like the well studied HD 209458b (ISwain et al. 



Snellen et al. 



20101 ). 



200s 



The properties (composition, temperature/pressure profiles, etc.) of exoplanet 
atmospheres retrieved from the emergent flux or from transit spectroscopy depend on the 
molecular opacities used to simulate the observed spectra. Spectroscopic studies of brown 
dwarfs and exoplanet atmospheres therefore need improved line parameters for hot NH3 
and CH4. 

Currently the best line list for NH3 is the HITRAN 2008 database, which is incomplete 
and has insufficient hot band information for the astronomical applications considered here. 
Moreover recent theoretical work has d emonstrated that H I TRAN 2008 has many errors in 



quantum number assignments of NH3 (IHuang et al. 



2011b 



Yurchenko et al. 



20091). NH3 



ha s already been shown to contribute significantly to the SED of brown dwarfs at ~620 



K (jPelorme et al 



20081 ). Improved theoretical predictions for the vibration-rotation lines 
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will help to assign the complicated spectra of NH3 (IZobov et al 



2OIII ) presented in this 



paper. We provide 



ligh 



temperature line lists for NH3 using a similar approach to that of 



Nassar fc BernathI (120031 ): in addition we derive empirical lower state energies for many of 
the lines. The line lists can be used directly to model the SEDs of brown dwarfs leading to a 
better understanding of the T- /Y-dwarf boundary and perhaps allow the first identification 



of NH3 in an exoplanet atmosphere. 
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Hot NHj line list at 1000°C 




1400 1500 1600 1700 1800 1900 2000 2100 




1400 1500 1600 1700 1800 1900 2000 2100 

Wavenumber (cm ') 



Fig. 1. — A line spectrum comparison of NH3 between the HITRAN line list converted to 
1000°C (top), the cahbrated emission line hst at 1000°C (middle) and the final fine list 
containing 'hot' emission lines supplemented with HITRAN lines (with an intensity greater 
than 1 X 10-22 cm/molecule) at 1000°C (bottom). 
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EXPERIMENTAL METHOD 



La boratory emission s p ectra of hot NH3 were recorded using a system similar to that 



used by lNassar fc BernathI (|2003| ) to study hot CH4. A diagram of the experiment is shown 
in Figure [21 An alumina (AI2O3) sample tube is sealed with windows and evacuated. 
A constant slow flow of NH3 gas is introduced to the system and a constant pressure is 
maintained using a needle valve which helps to minimize the build up of impurities and 
loss of sample within the system. Surrounding the central 51 cm of the AI2O3 tube (121 
cm long) is a controllable tube furnace capable of maintaining stable temperatures up to 
1370°C. The radiation exiting the AI2O3 tube is focussed into a Fourier transform infrared 
(FT-IR) spectrometer using a lens and the distance between the window on the end of the 
AI2O3 tube and the FT-IR spectrometer entrance aperture was 20 cm. This volume was 
purged with 'dry air' (H2O absent) to minimize the presence of H2O absorption lines in the 
emission spectra. 




Fig. 2. — A diagram of the experimental setup used (not to scale). The FT-IR spectrometer 
is on the left and the AI2O3 tube is on the right, the dashed line indicates the line of sight of 
the FT-IR spectrometer and the labels refer to: A. Scanning corner-cube; B. Beamsplitter; 
C. Fixed corner-cube; D. Dry air purge; E. Lens; F. Window; G. NH3 exits to vacuum pump; 
H. Tube furnace; I. Stainless steel end cap with rubber o-ring; J. NH3 enters from cylinder; 
K. and L. Cooling H2O coils; M. Stainless steel end cap with rubber o-ring; N. Window; O. 
Dry air exit; and P. MCT detector. 

We recorded all spectra in two sections to improve the signal-to-noise ratio (see 
Figure [3]) and the experimental conditions for each region are listed in Table [21 The two 
selected regions were a consequence of the available filters and constituent materials of the 
system. The first region covered the range 740 - 1690 cm~^ and was recorded with thallium 
bromoiodide (KRS-5) windows, a potassium bromide (KBr) beamsplitter, a zinc selenide 
(ZnSe) lens and a liquid nitrogen cooled mercury-cadmium telluride (MCT) detector. The 
lower wavenumber cut off was limited due to the band gap of the MCT detector and a small 
part of the NH3 V2 umbrella mode could not be observed. 




Fig. 3. — The two observed infrared emission spectral regions of NH3 at 1000°C. 
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Table 2. Experimental conditions. 



Parameter 


Region 1 


Region 2 


Spectral Region (cm~^) 


740 - 1690 


1080 - 2200 


Detector 


MCT 


MCT 


Beamsplitter 


KBr 


KBr 


Windows 


KRS-5 


CaF2 


Lens 


ZnSe 


CaF2 


Scans 


240 


80 


Resolution (cm~^) 


0.01 


0.01 


Aperture (mm) 


3.15 


2.5 


NH3 pressure (Torr) 


5.0 


1.0 


ZerofiUing factor 


xl6 


xl6 
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The second region covered the range 1080 - 2200 cm~^ (ahhough no hnes were 
observed above 2100 cm~^) and was recorded with calcium fluoride (CaF2) windows, a KBr 
beamsphtter, a CaF2 lens and liquid nitrogen cooled MCT detector. This second region 
covered the NH3 1/4 bending mode, plus the R-branch of the z/2 mode. For both regions, 11 
spectra were recorded at 100°C intervals between 300°C and 1300°C. The twelfth spectrum 
was recorded at 1370°C because the furnace had difficultly reaching 1400°C. The resolution 
was chosen based on the pressure and Doppler broadening widths and the final line list was 
created by cutting both regions at 1350 cm~^ and combining them to make a list of lines 
ranging from 740-2100 cm~^ for each temperature. 

A thermocouple within the interior of the tube furnace combined with a programmable 
controller was used to measure the furnace temperature and maintain high temperatures 
over long time periods. It was necessary to water cool the ends of the AI2O3 tube to avoid 
damaging the rubber o-rings at both ends of the tube which maintained the pressure seal. 
As a result, the ends of the AI2O3 tube are at a lower temperature than the central portion 
and because of the strong temperature dependence, the observed infrared emission spectra 
are dominated by the high temperature central portion of the heated sample. 

All spectra contain emission lines of NH3 assumed to be at the furnace temperature. 
The baseline of each spectrum depends upon the temperature and arrangement of the 
system; it originates primarily from thermal emission from the AI2O3 tube walls. The 
spectra also contain a small number of absorption lines due to H2O vapor (see Figure [3]). 
We also observed absorption of strong lines due to cold NH3 in the beam path, mainly in 
the end of the tube nearest the spectrometer. A description of how we deal with these 
absorption lines is given below. 



Once the spectra had been recorded, the peak-picking programme WSpectra (ICarleer 



2OOII ) was used to identify all of the emission peaks at each temperature; the number 
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of lines in each file is given in Table [31 The same program was then used to fit a Voigt 
profile to every emission line to create a line list (wavenumber, arbitrary intensity) for each 
temperature. The lines in the overlapping 1080 - 1690 cm^^ portion were compared at 
each temperature and calibrated to Region 1. It was decided that the two regions would be 
joined at 1350 cm~^ as there is a natural minimum in the line emissions and this value was 
approximately in the center of the overlapping section. Hence from here onwards, Region 
1 contains lines in the range 740 - 1350 cm~^ and Region 2 contains the lines 1350 - 2100 
cm~^. 
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Table 3. The number of measured emission lines. 



Temperature (°C) 


Region P 


Region 2^ 


300 


2,924 


4,014 


Ann 






500 


5,796 


7,186 


600 


7,665 


7,433 


700 


8,922 


9,810 


800 


8,841 


8,562 


900 


10,090 


10,851 


1000 


12,712 


9,526 


1100 


12,665 


9,752 


1200 


11,174 


9,992 


1300 


8,929 


9,233 


1370 


7,405 


6,366 


'^The upper (lower 
was chopped at 1350 


) limit of 
cm~^. 


region 1 (2) 
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A system response function for both regions was needed to account for the contribution 
from the system setup (e.g. windows, beamsphtter, lens, filter and detector) which alters 
the strength of the lines observed. This was achieved by recording the blackbody spectrum 
emitted from a solid graphite rod placed at the center of the AI2O3 tube; the end facing the 
spectrometer was machined into a concave cone to allow only blackbody emissions to be 
detected. The spectrum was then compared with a theoretical blackbody spectrum at the 
same temperature and normalized. The detector response curve for Region 2 and the effect 
on the line intensities can be seen in Figure HI 



- 19 - 




1200 1400 1600 1800 2000 2200 



Wavenumber (cm"^) 




1200 1400 1600 1800 2000 2200 

Wavenumber (cm'^) 



Fig. 4. — The top panel shows an observed uncahbrated NH3 hne spectrum at 1000 °C and 
the system response (blue) for region 2. The red dashed line indicates the low wavenumber 
cut-off for this region (i.e. 1350 cm~^). The bottom panel shows the resulting lines corrected 
for system response. 



All 24 fitted line lists (12 for each region) were wavenumber calibrated by comparison 
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20091 ). Strong, clear 



to 20 clean lines in the 2008 HITRAN database (iRothman et al. 
symmetric lines were favored and a calibration factor was then applied to correct the 
wavenumber scale. A typical calibration factor (obtained by dividing the HITRAN line 
center by observation) for region 1 was 1.000001739 (from the 900 °C line list) resulting 
in a shift of 0.00209 cm~^ at 1200 cm~^ and a typical calibration factor for region 2 was 
1.000000990 (from the 900 °C line list) resuhing in a shift of 0.00158 cm~^ at 1600 cm'^ 
The overall accuracy of our wavenumber scale is better than ±0.002 cm~^ after calibration. 



THEORY 



Our analysis method is not unique and has successfully been applied b y many workers 



most recently in the sub-millimeter regime to study astrophysical 'weeds' ( iFortman et al. 



20101). It involves a comparison of the observed intensities to a reliable data set (in our 
case, the 2008 HITRAN database). In order to make our e mission line lists compa. rable 



to HITRAN we converted them into absorption intensities (INassar fc Bernath 



20031 ). The 



relationship between emission and absorption intensities is given by 



5-, 



S, 



emission 



absorption 



where S^mission is the intensity of the emitted line, u is the frequency and T is the 
temperature. 

The HITRAN database provides line intensities at 296 K so they were convert ed to the 



relevant temperatures of the fitted line lists. The intensity of a line (S") is defined (IBernath 
20051 ) in SI units as 



, 27i^uioSj'j>' 
o = — z — ; — - — exp 



E 



Low 



1 — exp 



hu 



10 



kT 



(2) 



3eohcQ \ kT 

where Sj'j» is the square of the transition dipole moment, Q is the total internal partition 
function and E^ow is the lower state energy of the line. If the only changing variable is the 
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temperature, then the equation can be used to obtain the intensity of the same hne at any 
new temperature by dividing by a reference intensity (S'q), giving 



S' Qo 
— = — exp 



5, 



Q 



Et 

-'-'Low 



E 



Low 



kT 



l-exp(-^) 



exp 



huio 
' kTo 



(3) 



Sq and Qo refer to the reference temperature Tq, and S' and Q are at the new temperature 
T. The total internal partition functions for temperatures between 70 - 3005 K can be 
calculated by using the empirical equation 



Q{T) = a + bT + cT^ + dT^ 



(4) 



as given by 



Gamache et al. 



( 120001 ) ■ The constants in this equation change depending on the 



temperature range and the calculated internal partition function for NH3 is given in Table 



SI 
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Table 4. The partition function of NH3. 



T(K) 


Partition Function^ 


296 


1,729.30 


573 


5,305.19 


673 


7,218.42 


111 


o fine on 


873 


12,665.85 


973 


16,582.17 


1073 


21,548.31 


1173 


27,755.32 


1273 


35,394.26 


1373 


44,656.18 


1473 


55,732.16 


1573 


70,299.76 


1643 


81,868.17 



'^Calculated using Equation 
4 tak en from iGamache et al. 
(I2OO0I). 
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The fitted line lists were compared with the transformed HITRAN line lists at the 
corresponding temperature. This allowed an intensity linear calibration factor to be 
determined from a comparison of our intensities to HITRAN. These calibration factors 
are summarized in Table [5] and result in 24 intensity calibrated line lists (one for each 
temperature per region). The two regions were then combined to provide 12 line lists at 
each temperature. 

All calibrated line lists were then compared (to within ±0.005 cm~^) to identify the 
same lines at every temperature. Rearranging Equation 3 and taking the natural log of 
both sides gives 



In 



SQ 
SqQo 



1 — exp 



kTo I 



l_exp(-^) 



f SQRq 
\SoQo 



R 



E 



Low 



E 



Low 



kTo 



kT 



(5) 



so a plot of \n{SQRo/ SqQqR) against 1/kT yields a straight line. The slope of the 
fitted line provides the empirical lower state energy (E^ow) as illustrated with five typical 
lines which are present in all 12 line lists (Figure [5]). Since all of the lines are not present in 
HITRAN, only the Elow values of z>(4) [Elo«,(4) = 1564.0 cm'^] and i>(5) [^lo«,(5) = 1808.6 
cm^^] can be compared to HITRAN which gives an accuracy of 7% and 2% respectively 
(compared to HITRAN E^ow values of 1673.5 and 1771.3 cm~^). Note that we have chosen 
to work with the same units as in the HITRAN database. Conversions to other units, 



particularly for intensities (e.g. osci 
obtained by using equations given in 



lator str e ngths , / and Einstein A values) can be 



Bernathl (|2005|) 
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Table 5. Linear calibration factors (/) to convert the measured line intensities (Sarb) into 

HITRAN units, S' in (cm/molecule). 



Temperature (°C) 


Region l'^ 
/i (xl0i9) 


Region 2*^ 
/2 (xlO^^) 


300 


1.0 


15.0 


400 


1.0 


10.0 


500 


1.0 


7.0 


600 


1.0 


5.0 


700 


1.25 


4.0 


800 


1.25 


3.5 


900 


1.25 


3.0 


1000 


1.4 


2.3 


1100 


1.5 


2.8 


1200 


3.5 


5.5 


1300 


4.0 


7.0 


1370 


4.5 


7.0 



^Where S' = fi^2)Sarb- 
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Fig. 5. — The natural log of the intensity ratios of five lines present in all 12 line lists are 
plotted against 1/kT (in cm~^), where T is the temeprature. The slope of each line directly 
yields the empirical lower state energy and for the lines plotted, Elow{^) = 2519.7 cm~^ 
[z>(l)], Elo^{2) = 2091.7 cm-i [i>(2)], Elow{3) = 2221.7 cm~^ [i>(3)], Elo»(4) = 1564.0 cm"! 
[z>(4)], and ^lo«;(5) = 1808.6 cm"! [z>(5)]. Points above 1000°C {1/kT < 0.00113 cm"!) were 
neglected in the fit. 
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4. RESULTS AND ANALYSIS 

We were able to obtain the lower state energies for the majority of the observed lines 
and a plot of the results in the region 740 - 2100 cm~^ can be seen in Figure [61 There is 
good agreement with the HITRAN lower state energies for the same region displayed in 
Figure [TJ and the pattern of lines is similar. Note that because of self-absorption cold NH3 
in the end of the AI2O3 tube means that the strong room temperature lines are missing. 
They have been added back into the final line lists from HITRAN to make each line list 
complete. 



5000 



4000 
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2000 



1000 




800 



1000 



1200 



1400 



1600 



1800 



2000 



Wavenumber [cm" ] 



Fig. 6. — Empirical lower state energies [Elow) determined from all line lists. The colour 
refers to the number of temperatures at which the line appeared (12 being the maximum). 
It is indicative of the quality of the calculated E^ow 
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Fig. 7. — Lower state energies taken from HITRAN. P- Q- and R-branches can clearly be 
identified in the plot. 

The majority of the i>2 bending mode (umbrella mode) can be seen in Figures [6] and [71 
The band gap of the MCT detector means that we cannot detect any lines below 740 cm~^. 
What is most striking about this region is the strong Q-branches rising sharply between 900 
and 1000 cm~^. Our calculated empirical lower state energies only identify these Q-branches 
above Elow ~ 1500 cm~^. The patterns in Figures [6] and [7] can be explained qualitatively 
using the rigid rotor energy level formula, 



E.j^K = E, + BJ{J + 1) + {C- B)K^ 
E^,j,K[cm-i] = E^ + 9.94J( J + 1) - 3.75i\:^ 



(6) 
(7) 
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with the vibrational energy. The parabohc features increasing to the right (left) of 
the Q-branches are due to the R-branches (P-branches) of different vibrational levels. The 
small nearly vertical strips within the 1^2 region around 950 cm~^ are due to K-structure 
(A J = ±1, AK = 0) and the doubling of the band due to inversion of the NH3 molecule is 
clearly evident. Also apparent at high Elow in Figure [7] are decreasing parabolic features 
due to forbidden K-structure {AK = ±3, AK = ±6) transitions. The R-branch in the z/2 
region at Elow ~ 1800 cm~^ above the ground state is likely due to the population in the 
f 2 = 2 vibrational level (i.e., 3z/2 — 2z/2 band). 

The Q-branch of the 1^4 bending mode can also be seen rising sharply at 1650 cm^^; 
the P- and R-branches are also visible as gradually increasing parabolas from the same 
point. The majority of lines in the HITRAN database are fundamental transitions but the 
1^2 + — 1^2 hot band is also present at 1600 cm^^ above an Elow ~ 1000 cm^^. Our high 
temperature observations include several additional hot band transitions and a hot band is 
clearly visible in the z/4 region (Figure [6]) with an Eiow ~ 1500 cm^^ which can likely be 
assigned to the 2z/4 — 2/4 hot band. 

What is clearly illustrated in Figures [6] and [71 is the lack of observed transitions with 
lower state energies below approximately 1500 cm"^ for the z/2 band and ~ 1000 cm^^ for 
the z/4 band. The main reason for this is due to the experimental setup (shown in Figure 
[2]). The NH3 is heated at the center of the tube and emission occurs in this hot region; this 
emitted radiation travels along the tube to the window (maintained at room temperature). 
NH3 close to the window is at approximately room temperature and absorbs some of the 
emitted radiation, mainly for the low J fundamental lines. When the emission is observed, 
the majority of emission lines are from hot bands as well as high J fundamental lines that 
are not populated at room temperature. 

To make our line lists complete we have inserted all HITRAN lines with an intensity 
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above 1 x 10"^^ cm/molecule at each temperature. As some of the observed emission 
hnes are due to these HITRAN hues, the observed hues have been removed to avoid the 
same hne appearing twice. Table E] hsts the number of hnes in each file broken down into 
added HITRAN lines and observed lines (the total number of observed lines starts to 
decrease above 1100 °C because of NH3 decomposition in the cell). The consequence of this 
procedure can be seen from a lower state energy plot of the 1000°C line list (Figure [H]). The 
12 line lists at each temperature have been amalgamated into one table which is available 
in a machine readable version online. A sample of this final line list for all temperatures is 
provided in Table [71 
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Table 6. The number of lines present in each Une Ust. 



Temperature (°C) 


Total Lines 


Observed Lines 


Added HITRAN Lines 


300 


8,102 


3,611 


4,491 


400 


10,264 


3,758 


6,506 


500 


14,043 


3,816 


10,227 


600 


16,305 


3,784 


12,521 


700 


19,965 


3,688 


16,277 


800 


19,003 


3,574 


15,429 


900 


22,680 


3,431 


19,249 


1000 


24,168 


3,295 


20,873 


1100 


24,506 


3,173 


21,333 


1200 


23,151 


3,084 


20,067 


1300 


20,089 


2,955 


17,134 


1370 


15.410 


2.873 


12.537 
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Fig. 8. — Lower state energies plotted from the 1000°C line list. The colour refers to the 
quality of the calculated empirical lower state energy. 

90 



Table 7. A sample of 10 lines from the 1000°C line list. 



Temperature (°C) Wavenumber (cm ^) Intensity (cm/molecule) Elow (cm ^) Quality Identifier 



1000 


1591.8243452 


0.68853731E- 


■22 


0.1 584654 lE+04 


2 


1000 


1591.8273898 


0.68853731E- 


•22 


O.OOOOOOOOE+00 





1000 


1591.8536200 


0.62000000E- 


■21 


0.14309118E+04 


IH 


1000 


1591.9071469 


0.36946435E- 


■21 


0.26669918E+04 


2 


1000 


1591.9233494 


0.42613992E- 


■21 


0.84666645E+03 


3 


1000 


1591.9558781 


0.41360202E- 


■21 


0.15639455E+04 


1 


1000 


1591.9876537 


0.16819148E- 


•21 


0.33495921E+04 


2 


1000 


1592.0143421 


0.75043627E- 


■22 


0.30381832E+04 


3 


1000 


1592.0369615 


0.79134180E- 


■22 


0.15844298E+04 


2 


1000 


1592.0767117 


0.97189882E- 


■21 


0.16858822E+04 


1 
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The estimation of the error in the empirical lower state energies is difficult. As a 
measure of the accuracy of our empirically-determined lower state energies, we determined 
the differences between our values and those given in HITRAN for the lines that are present 
in both data sets. The average Eiow difference was found to be 147 cm~^ which corresponds 
to 3.5%. However, the average magnitude of the deviation is determined to be 364 cm~^ 
corresponding to 30%. Since these lines have been largely removed from the final line lists 
and replaced by the actual HITRAN values due to self absorption issues, we consider these 
errors to be an over estimation. 

Another measure of the quality of the empirical lower state energy can be obtained 
from the number of points (i.e. the number of line lists each particular line appears in) used 
to calculate the value. For lines occurring in 10 - 12 spectra the Elow was calculated from 
the points below 1000°C and has a quality factor of '1', these are considered to have the 
least error. For lines in which an Elow can still be determined from points below 1000°C 
(i.e. occur in at least 3 spectra below 1000°C) but occur in less than 10 spectra overall, 
the quality factor is '2'. For any line in which an Elow is still able to be determined but 
only from 3 points including spectra above 1000°C, then the quality factor is '3', and are 
considered to have the greatest error. Any line where Elow cannot be accurately determined 
(i.e. less than 3 points), then the quality factor is '0'. The HITR AN lines which o ccur in 
the line lists have a quality factor of 'IH' and are are taken from iRothman et al.l (120091 ). 
The quality factors can be seen in Table [7] and the effect on the lower state energy can be 
seen in Figure [HI 



5. DISCUSSION 



The ammonia line lists have been wavenumber calibrated with the HITRAN 2008 
database and are accurate to within ±0.002 cm~^. The intensity was also calibrated 
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with the HITRAN 2008 database by extrapolating the intensities using Equation 3 to the 
temperatures of the recorded spectra. An error in the cahbrated intensity of perhaps a 



factor of two is suggested for our hn e hsts. Although t 



are notoriously difficult to calibrate ( iNassar fc Bernath 



lis see ms large, emission intensities 



20031). 



We have accounted for the instrument response by recording a blackbody spectrum 
and correcting the intensities with the instrument response function. The most notable 
corrections were obvious towards the edge of each spectral region (see Figure H]). Blending 
of lines also causes errors in line intensities and in the empirical lower state energy. These 
blended lines cause problems when matching our observations with HITRAN and make it 
difficult to identify the contribution of each line to the observed peak. We have tried to 
reduce our intensity calibration error by comparing as many spectral lines with HITRAN 
lines as possible. Each comparison gives an intensity calibration factor and after outliers 
have been removed a general factor can be obtained. The temperature of the experiment 
was accurate to within ±10 K but this only applies to the central portion of the AI2O3 
tube which is surrounded by the furnace. There is therefore a temper ature gradient with in 



the A I9O3 tube which affects the shape of some lines as discussed by 



Nassar fc Bernath 



( I2OO3I ). This effect can be observed in the emission spectra as we see absorption of the low 
J fundamental NH3 lines due to cooler gas towards the tube ends change to emission as the 
temperature of the tube increases (see Figure [9]). 
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Fig. 9. — A short (1 cm~^) section of the NH3 emission spectra. Notice that the hne at 
1095.13 cm~^ (1^(7) of the fundamental V2 band) changes from absorption to emission as the 
temperature increases. 

At higher temperatures the AI2O3 tube ends are sufficiently hot and the emission from 
the center is stronger so they no longer completely absorb the low J fundamental lines 
leading to differences between successive spectra starting around 1000 °C. As long as the 
center of the tube is 'cool' (<900 °C) the fundamental low J lines appear in absorption, 
however for temperatures greater than 1000 °C the fundamental low J lines start to appear 
in emission. Moreover, because the Doppler broadening depends on temperature the 
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emitted lines are broader than the absorption hnes. This leads to an absorption feature in 
the middle of an emission line which can sometimes cause apparent splitting of the line 
(e.g. the line at 1095.13 cm~^ in Figure [H]). To solve this problem, we have processed each 
line list and inserted HITRAN lines with intensities greater than 1 x 10~^^ cm/molecule 
at that temperature, denoted by 'IH' in the line lists. We also simultaneously removed 
the emission lines showing this splitting effect, and also any lines which are within ±0.005 
cm~^ of each HITRAN line with an intensity greater than 1 x 10~^^ cm/molecule. Our 
line lists therefore provide a complete representation of NH3 between 740 - 2100 cm~^ at 
temperatures between 300 - 1370°C. 

A useful check for the quality of the calibrated intensities determined between 740 - 
2100 cm~^ can be ve rified by an intensity su m comparison with the temperature scaled 



HITRAN intensities ( Nassar fc Bernath 



20031 ). At 300°C, the sum of the intensities in our 



linelist is only 7% greater than the sum of the intensity scaled HITRAN lines but as the 
temperature increases we continue to exceed the HITRAN total intensity with a maximum 
of 105% greater at 1200°C (Table E]). We attribute our larger values to an increase in 
the number of hot band lines present in the spectrum. These hot lines only have a small 
contribution at 'low' temperatures (300°C) but are more significant for higher temperatures 
(> 1000°C) as hot transitions become more favorable. Indeed, this comes as no surprise 
because the HITRAN database consists mainly of fundamental transitions and highlights 
how our line lists are more suitable for astrophysical applications. 
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Table 8. Intensity sums for each line list and HITRAN between 740-2100 cm~^ 



Temperature Intensity Sum^ HITRAN intensity sum Ratio to 
(°C) (cm-i) (cm-i) HITRAN 



"inn 


O 

L. 


. 1 u 


X 




-17 
















300 


2. 


.24 


X 


10- 


-17 


2, 


.39 


X 


10- 


-17 


1, 


.07 


400 


2. 


.04 


X 


10- 


-17 


2, 


.27 


X 


10- 


-17 


1, 


.11 


500 


1, 


.82 


X 


10- 


-17 


2 


.12 


X 


10- 


-17 


1, 


.16 


600 


1, 


.59 


X 


10- 


-17 


1 


.95 


X 


10- 


-17 


1, 


.23 


700 


1, 


.37 


X 


10- 


-17 


1 


.85 


X 


10- 


-17 


1, 


.35 


800 


1, 


.17 


X 


10- 


-17 


1 


.68 


X 


10- 


-17 


1, 


.44 


900 


9, 


.89 


X 


10- 


-18 


1 


.54 


X 


10- 


-17 


1, 


.56 


1000 


8, 


.34 


X 


10- 


-18 


1 


.49 


X 


10- 


-17 


1, 


.77 


1100 


7, 


.03 


X 


10- 


-18 


1 


.33 


X 


10- 


-17 


1, 


.89 


1200 


5, 


.94 


X 


10- 


-18 


1 


.22 


X 


10" 


-17 


2, 


.05 


1300 


4, 


.92 


X 


10- 


-18 


9 


.64 


X 


10- 


-18 


1, 


.96 


1370 


4, 


.35 


X 


10- 


-18 


8, 


.05 


X 


10- 


-18 


1, 


.85 



*The intensity sums were calculated before the removal and addition of 
lines from the calibrated line lists. 
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A similar metho d to determine empirical lower state energies has been used by 



Fortman et al 



(120101 ) by recording numerous millimeter wave spectra with the Fast Scan 
Submillimeter Spectroscopic Technique (FASST) as a f unction of temperat u re. F or NH3, 



Brown fc Margolid (119961 ) for 



empirical lower state energies were also determined by 
unassigned lines in the 4791 - 5294 cm~^ region. Another method has been used recently 
to o btain lower state e nergies for near infrared lines of CH4 by cooling an absorption 



cell (IWang et al 



20101). For the FT-IR method presented here, we have line lists at 12 
individual temperatures each of which requires many hours to obtain. As a result, we are 
only able to use a maximum of 12 points in the calculation of the empirical lower state 
energy for each line and in fact, this was reduced to 8 points as the intensity of the lines 
above 1000°C begin to decrease in strength due to the dissociation of the molecule. Our 
accuracy depends on the number of points used in the calculation of Elow and for this 
reason, our line lists also include a quality factor which only applies to the empirical lower 
state energy. This Elow method is used to determine an important characteristic of each 
line whilst knowing nothing of the spectroscopic quantum number assignments. 

The inclusion of HITRAN lines in our line lists has benefits as well as problems. The 
main benefit is that the line lists now contain the strong lines which were absent from 
the observations and as a result are more suitable for simulation of astronomical spectra 
(Figured]). As a consequence of including these strong lines (1 x 10^^^ cm/molecule), some 
observed lines have been removed because they have the same frequency as the HITRAN 
line. It was decided that the benefits of including the HITRAN lines outweighed the costs 
of removing coincident emission lines. 

The errors associated with the empirical lower state energies have been based on 
comparisons with HITRAN however, as the majority of lines used to determine the errors 
are self absorbed, it does not accurately represent all of the lines in each spectrum. Indeed, 
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at low temperatures the strong lines which are present display self absorption even if it is 
not immediately apparent from the line shape and at high temperatures there are few lines 
to draw conclusive comparisons from, since they are not given in HITRAN. The overall 
maximum error of the empirical lower state energies in the line list should therefore be taken 
as 30%, but we believe that the quality factor gives a better indication of how accurate the 
empirical lower state energy is for each line. 

There are two approaches to obtaining an ammonia line list to model brown dwarf 
and exoplanet spectra: the first is experimental measurement consistent with the work 
presented here and the second is an ab initio theoretical calculation. These two methods 
each have their own strengths and weaknesses, but in general are complementary. The 
experimental approach provides fewer lines but more accurate line positions, while the 
theoretical approach provides many more lines but with reduced line position accuracy. 
Calculated line lists are particularly useful for assigning quantum numbers to the measured 
lines, which can then be used to improve the calculations. Such l ine assignments o f new 



NH3 data are currently underway and will be reported elsewhere (jZobov et al.ll201ll ). 



The line lists provided are given in the format shown in Table [7] and can be obtained 
as electronic supplements with extra details and spectra provided upon request. It is 
important to note that the line lists are temperature specific and should only be used at the 
appropriate temperature. If the intensities need to be scaled, then the closest temperature 
line list should be chosen and the empirical lower state energies can be used to adjust the 
lines intensities as needed (using Equation 3) to reach additional temperatures. 

We are currently in the process of completing our NH3 measurements in the 3 micron 
region and are set to record new spectra in the near infrared; a similar comprehensive study 
of CH4 is also underway. 
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6. CONCLUSION 

We have produced 'hot' experimental hne hsts of NH3 at temperatures 300 - 1300°C in 
100°C intervals and at 1370°C. Each hne hst has been wavenumber and intensity cahbrated 
using the HITRAN 2008 database and from this we have been able to derive empirical 
lower state energies from the temperature dependance of line intensities. The line lists can 
be used directly to model NH3 between 740 - 2100 cm~^ (4.76 - 13.5 iim) in brown dwarfs 
and exoplanets. 

Support for this work was provided by a Research Project Grant from the Leverhulme 
Trust and a Department of Chemistry (University of York) studentship. 
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